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►-_ ' Abstract. Precession and nutation of the Earth depend on the Earth's dynamical flattening, H, which is closely 

^H , related to the second degree zonal coefficient, J2 of the geopotential. A small secular decrease as well as seasonal 

QQ ■ variations of this coefficient have been detected by precise measurements of artificial satellites (Nerem et al. 

^SJ ' 1993, Cazenave et al. 1995) which have to be taken into account for modelling precession and nutation at a 

microarcsecond accuracy in order to be in agreement with the accuracy of current VLBI determinations of the 
Earth orientation parameters. However, the large uncertainties in the theoretical models for these J2 variations (for 
example a recent change in the observed secular trend) is one of the most important causes of why the accuracy 
of the precession- nutation models is limited (Williams 1994, Capitaine et al. 2003). We have investigated in this 
O ' paper how the use of the variations of J2 observed by space geodetic techniques can influence the theoretical 

expressions for precession and nutation. We have used time series of J2 obtained by the "Groupe de Recherches 






H 



^ ■ en Geodesic spatiale" (GRGS) from the precise orbit determination of several artificial satellites from 1985 to 2002 

' to evaluate the effect of the corresponding constant, secular and periodic parts of H and we have discussed the 

best way of taking the observed variations into account. We have concluded that, although a realistic estimation 

of the J2 rate must rely not only on space geodetic observations over a limited period but also on other kinds of 

^ ' observations, the monitoring of periodic variations in J2 could be used for predicting the effects on the periodic 

•/j , part of the precession-nutation motion. 

. ' I Key words, astrometry - reference systems - ephemerides - celestial mechanics - standards 

T-H ■ 

I ■ 1. Introduction cent models for precession (Williams 1994, Capitaine et 

(-^ ' al. 2003) or nutation (Souchay & Folgueira 1999, Mathews 
Vi : Expressions for the precession of the equator rely on ^^ ^^ ^^^^^ Lambert & Capitaine 2004) in which either 
values for the precession rate in longitude that have ^^^ ^^^^j^^ ^^ ^^^ ^^^^1 variations of this coefficient are 
been derived from astronomical observations (i.e. obser- ^^pii^jtly considered through simplified models, 
vations that were based upon optical astrometry until the 
C3 ; lAU 1976 precession, and then on Very Long Baseline The recent implementation of the lAU 2000 A 
Interferometry (VLBI) observations for more recent mod- precession-nutation model guarantees an accuracy of 
els). The lAU 2000 precession-nutation model provided about 200 ^as in the nutation angles, and all the pre- 
by Mathews et al. (2002) (denoted MHB 2000 in the dictable effects that have amplitudes of the order of 10 ^as 
following), that was adopted by the lAU beginning on have therefore to be considered. One of these effects is the 
1 January 2003, includes a new nutation series for a non- influence of the variations (AH) in the Earth's dynami- 
rigid Earth and corrections to the precession rates in Ion- cal flattening, which are not explicitly considered in the 
gitude and obliquity that were estimated from VLBI ob- lAU 2000A precession-nutation model. Furthermore, the 
servations during a 20-year period. The precession in Ion- lAU 2000 precession is based on an improvement of the 
gitude for the equator being a function of the Earth's dy- precession rates values derived from recent VLBI measure- 
namical flattening H, observed values for this precession ments, but it does not improve the higher degree terms in 
quantity are classically used to derive a realistic value for the polynomials for the precession angles ipA, (^a of the 
H. Such a global dynamical parameter of the Earth is equator (see Fig.p. This precession model is not dynami- 
generally considered as a constant, except in a few re- cally consistent because the higher degree precession terms 
are actually dependent on the precession rates (Capitaine 
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solutions due to the limited span of the available data 
(Capitaine et al., 2004). One alternative way for such an 
improvement is to improve the model for the geophysical 
contributions to the precession angles and especially the 
influence of AiJ (or equivalently AJ2). 

The H parameter is linked to the dynamical form- 
factor, J2 for the Earth (i.e. the C20 harmonic coefficient 
of the geopotential) which is determined by space geodetic 
techniques on a regular basis. Owing to the accuracy now 
reached by these techniques, the temporal variation of a 
few Earth gravity field coefficients, especially AC20, can 
be determined (for early studies, see for example Nerem et 
al. (1993), Cazenave et al. (1995) or Bianco et al. (1998)). 
They are due to Earth oceanic and solid tides, as well 
as mass displacements of geophysical reservoirs and post- 
glacial rebound for AC2o. This coefficient C20 can be 
related to the Earth's orientation parameters and more 
particularly to the Earth precession-nutation, through H . 
The purpose of this paper is to use space geodetic deter- 
mination of the geopotential to estimate Aif , in order to 
investigate its influence on the precession-nutation model. 
The C20 data used in this study have been obtained from 
the positioning of several satellites between 1985 and 2002. 
We estimate also the constant part of H , based on such 
space geodetic measurements, and compare its value and 
influence on precession results with respect to those based 
on VLBI determinations. 

In Sect. 2 we recall the equations expressing the equa- 
torial precession angles as a function of the dynamical flat- 
tening H. We provide the numerical values implemented 
in our model, compare the values obtained for H by var- 
ious studies and discuss the methods on which they rely. 
In Sect. 3 the relationship between Ai/ and AC20 is dis- 
cussed, depending on the method implemented. We ex- 
plain how these geodetic data are taken into account in 
Sect. 4. We present our results in Sect. 5, and discuss them 
in the last part. We investigate how the use of a geodetic 
determination of the variable geopotential can influence 
the precession- nutation results, considering first the pre- 
cession alone, and second the periodic contribution. 

In the whole study, the time scale for t is TT Julian 
centuries since J2000, which will be denoted cy. 

2. Theoretical effect of /^H on precession 

This section investigates the theoretical effect of the vari- 
ations AiJ in the Earth's dynamical flattening on the pre- 
cession expressions. 

2.1. Relationship between H and the precession of the 
equator 

The two basic angles V'A and uja (see Fig.[T|) for the preces- 
sion of the equator are provided by the following differen- 
tial equations (see Eq. (29) of WiUiams (1994) or Eq. (24) 
of Capitaine et al. (2003)): 

dil^A 



Mean equator of date 




Ecliptic of date 



Ecliptic of 
epoch J2000 



Mean equator of epoch 



Fig. 1. Angles ipA and uja for the precession of the equa- 
tor: 7m is the mean equinox of the date and 70 is the 
equinox of the epoch J2000.0. 



duJA 

dt 



re cosxA + ir^, sine a) sinxA 



(1) 



where r^ and r^ are respectively the precession rates in 
longitude and obliquity, eA is the obliquity of the eclip- 
tic of date and xa the planetary precession angle, de- 
termining the precession of the ecliptic. Updated expres- 
sions for these precession quantities are given in Capitaine 
et al. (2003). An expression for the precession rates, r^ 
in longitude and r^ in obliquity, is provided in detail in 
Williams (1994) and Capitaine et al. (2003) as a function 
of various contributions. The precession rate in longitude 
can be written as r^ = tq + ri t + r2 t'^ + r^ t^ where 
the largest first-order term in vq is the luni-solar contribu- 
tion denoted /oiks coseoi where eo is the obliquity of the 
ecliptic at J2000. It is such that (Kinoshita 1977, Dehant 
& Capitaine 1997): 



/oi|is = k,n Mo + kg So 



(2) 



in which Mq and 5o are the amplitudes of the zero- 
frequency Moon and Sun attractions, respectively, and: 

+ rriQ F2^ ^ 



3 H 



H K„ 



3 H 



m© 



H K, 



(3) 
(4) 



svauJA 



dt 



{r^ sine a) cosxA-n sinxA 



rriQ + mm + m-© ^ 

In the above expressions, H is the Earth's dynamical flat- 
tening, m„i, TTiQ and to^ are the masses of the Moon, the 
Earth and the Sun, respectively, rXm is the Moon mean 
motion around the Earth, uq the Earth mean motion 
around the Sun, fl is the mean angular velocity of the 
Earth and F2 a factor for the mean distance of the Moon. 
Current numerical values for such a problem are (Souchay 
& Kinoshita, 1996): 

Mo = 496303.66 x 10"^ 
5*0 = 500210.62 x 10"'^ 

fc„ = 7546". 7173289 /cy (5) 

fc, = 3475". 1883295 /cy 
/oi|^^coseo = 5040".6445/cy 
and (see Kinoshita, 1977): 
F2 = 0.999093142 
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Hence, the link between the precession of the equator 
{ipA and UJA angles) and the Earth's dynamical flattening 
(H) is shown by Eq. H]), Eq. ©, Eq. ©, Q and Eq. © 
of Sect. 2.2. Classically, H is related to /oiu c derived from 
observations by: 



H 



/oi 



\ls 



K^ Mo + Ks So 



(6) 



(7) 



2.2. Astronomical determination of H 

We can write rg as: 

ro = /oiiigCoseo + /oi|pi coseo 

+ H X lunisolar second order effects 

+ H X ( J2 and planetary) tilt effects 

+ J4 lunisolar effect 

— geodesic precession 

+ non- linear effects (Mathews et al., 2002) 



where /oiipj^ is the first order term of the planetary con- 
tribution (also proportional to H). Classically, H is de- 
rived from observationally determined values of rg. The 
measurement of rg should be corrected by removing the 
modelled contributions other than the lunisolar first order 
effect (see Eq. Q). Hence, we obtain a value for /oiUs7 
which is the only term with sufficiently large amplitude 
(of the order of 5000"/cy) to be sensitive to small changes 
in the value of the dynamical ellipticity H of the Earth 
(see Eq. (|6])). So, given the other contributions provided 
by the theory, we can derive the value of H from the ob- 
served value of To and the model for the lunisolar first 
order effects. 

A major problem consists in choosing the constant 
value of H. Indeed, depending on the authors, it differs 
by about 10~^ (Table[T]). This is due to the different mea- 
surements and models implemented (see Fig. 1 of Dehant 
& Capitaine (1997); Fig. 5 of Dehant et al. (1999)). On 
the one hand the optical measurements give values of the 
general precession in longitude pA referred to the eclip- 
tic of date, whereas VLBI gives measurements relative to 
space. On the other hand, the various constants and mod- 
els used for obtaining the value for H from a measured 
value (optical. Lunar laser ranging or VLBI) are different 
depending on the study considered (see Eq. (O). 

Classically, ipA is developed in a polynomial form of 
t as: ipA = -ipo + i'l t + i^2 f^ + 'fpa t^ ■ In Table [U we 
recall the different values used (i) for i/ji (i.e. the precession 
rate in longitude, ipi = tq), directly obtained from VLBI 
measurements, and (ii) for pi which is the observationally 
determined value of precession in the optical case: i/'i = 
Pi + Xi cos eo (Lieske et al., 1977). 

The computation of the lAU 2000 precession-nutation 
model by Mathews et al. (2002) is based on a new method 
which uses geophysical considerations. They adjust nine 
Basic Earth Parameters (BEP), including the Earth dy- 
namical flattening H . 



2.3. IVIethod and parameters used in this study 

Based on the paper by Capitaine et al. (2003), denoted 
hereafter P03, we use differential equations ([T]) in which 
H has been replaced hy H + AH (using Eq. (H]), Eq. ©, 
Eq. gl) and Eq. ^). We start from the P03 initial values 
for the variables uja, V'a, ca, Xa and pA, that are repre- 
sented as polynomials of time and rely on the numerical 
values given in Table[2] We solve Eq. ([T]) together with the 
other precession equations (e.g. see Eq. (26) and Eq. (28) 
of P03) with the software GREGOIRE (Chapront, 2003) 
that can process Fourier and Poisson expressions. We it- 
erate this process until we obtain a convergence of the 
solution. 



Table 2. Numerical values used in this study. H, -01 and 
LUi are integration constants. 



Initial values at J2000.0 


H 




^MHB = 3.27379492 x 10"^ 


Vi 




5038".481507/cy 


Wl 




-0".02575/cy 


Pi 




5028".796195/cy 


XI 




10".556403/cy 


eo 




84381".406 = 23°26'21".406 


Contributions 


to the 


precession rate in longitude (in "/cy) 



Lunisolar first order 5494.062986 x cos eo ^ 5040.7047 
Planetary first order 0.031 
Geodesic precession —1.919882 



3. Relationship between C20 and H 

3.1. Relation 

From the geodetic C20 variation series we can derive the 
corresponding variations of the dynamical flattening H. 
Indeed, knowing that J2 = — C20 = — \/5 C20, in the 
case of a rigid Earth, we can write (see Lambeck, 1988): 



H = C- 



A + B 



/C = 



M R, 



C 

M Re 

c~ 



J2 



(8) 



C2. 



M Re 

c 



C, 



20 



where A, B and C are the mean equatorial and polar mo- 
ments of inertia of the Earth. Al and Re are respectively 
the mass and the mean equatorial radius of the Earth. C20 
is the normalized Stokes coefficient (of degree 2 and order 
0) of the geopotential. 

But the Earth is elastic, so let us consider small varia- 
tions of 7f , C20 and the third principal moment of inertia 
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Table 1. Comparison between constants used for different determinations of the dynamical flattening {H): (1) the 
precession rate in longitude {ipi ) , (2) the speed of the general precession in longitude {pi ) , (3) the geodesic precession 
(pg) and (4) the obliquity of the ecliptic at J2000.0 (eo). The observational value actually used for each study is written 
in bold. 



Sources 


H 

(x 10^) 


(1) 

V'l 

( 


(2) 
Pi 

— in'Vcy 


(3) 
Pa 

) 


(4) 
eo 


Lieske et al., 1977 




5038.7784 


5029.0966 


-1.92 


23°26'21".448 


Kinoshita, 1977 and Seidelmann, 1982 


3.2739935 


5038.7784 


5029.0966 


-1.92 


23°26'21".448 


Williams, 1994 


3.2737634 


5038.456501 


5028.7700 


-1.9194 


23°26'21".409 


Souchay & Kinoshita, 1996 


3.2737548 


- 


5028.7700 


-1.9194 


23°26'21".448 


Bretagnon et al., 1997 


3.2737671 


5038.456488 


5028.7700 


-1.919883 


23°26'21".412 


Bretagnon et al., 2003 


- 


5038.478750 


5028.792262 


-1.919883 


23°26'21".40880 


Fukushima, 2003 


3.2737804 


5038.478143 


5028.7955 


-1.9196 


23°26'21".40955 


Capitaine et al., 2003 


3.27379448 


5038.481507 


5028.796195 


-1.919883 


23°26'21".406 


Mathews et al., 2002 


3.27379492 


5038.478750 


5028.7923 


-1.9198 


23°26'21".410 



of the Earth (C being its constant part and C33 its variable 
part). Then we obtain: 



H 



total 



M Re 

C 



1 



J2 



C33 ^2 total 



(9) 



we 



c-izjC being a small quantity of the order of 10 
consider the Taylor development of (1 4- c-^-ilCy . Then 
the total expression of H can be written as: 



H 



total 



MRe 

c 



J'. 



2 total 



1 



C33 
c 



m 



(10) 



where MR^^ /C x (C33/C)" J2 for n > 1 is smaller than 
10""'^"'^. So, in Eq. PH)) . considering (i) constant and vari- 
able parts separately and (ii) Eq. ([8]), we obtain: 



A„.^A., = _vi;i!flAft„ 



(11) 



where AJ2 = — AC20 = — V5 AC'20 corresponds to the 
variations of the Stokes coefficient J2. Generally, we can 
write: AJ2 oc C33/C (Lambeck, 1988). 

3.2. Computation of the ratio MR^ jC 

The coefficient MR^^ /C is usually obtained from the H 
and J2 values (see Eq. dS])). In order to determine the 
constant part of H , we can use (i) the Re, M and C values 
or (ii) the Clairaut theory (see Table [3]). 

First, recall the Earth geometrical flattening e: 



Re 



(12) 



where Rp and Re are respectively the polar and equatorial 
mean radii of the Earth. Second, recall the assumptions 
that the Earth (i) is in hydrostatic equilibrium and (ii) 
is considered as a revolutional ellipsoid. Hence, the first 



Clairaut equation gives the Earth geometrical flattening 
as a function of J2 and q. The approximations to the first 
and second order are respectively: 



:J2 



9^3 9 2 3 11 2 

J2 -\ J2 J2 Q Q 

2 2 8 14 ^ 56^ 



where the geodynamical constant is: 

LU^Re^ 



GM 



(13) 
(14) 

(15) 



= 3.461391 X 10"^ , lAG (Grotcn, 1999). 



Then, the following Radau equation can help us to deter- 
mine the expression of MRe^ /C: 



q/2 



H 



= i--yTTv = 



(16) 



where A is the d'Alembert parameter and 77 the Radau 
parameter, as: 



V 



5q 
2e 



(17) 



Hence, replacing e with Eq. ([13]) in Eq. (|T6l) and using 
Eq. ([5]) gives the Darwin-Radau relation as following: 



C 



MR? = 3-X = 3V-5'^^^ 



(18) 



Our tests have shown that Eq. p^ for e, in the expression 
p?)) of 77, gives more reliable results. 

In Table [3] we compare the various H values obtained. 
We denote (i) H* the value obtained with the Clairaut 
method and (ii) H** the value obtained using directly the 
Re, C and M values. Both are computed with Eq. ([8]) and 
a value for J2 of 1.0826358 x 10~^. Note that in contrast, 
lAG or MHB values (usually used) are determined from 
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astronomical precession observations and can be used to 
compute the C jMR,?' value. We can add that the differ- 
ences with Hmhb come from (i) the hydrostatic equilib- 
rium hypothesis in Clairaut's theory for the value H* and 
(ii) the poorly determined i?e, C and M values, for the 
value H** . This will introduce errors in the /S.H determi- 
nation, which we will study in Sect. 3.3. 

In the following, we will use the C / (^MRJ^) value de- 
termined with the Clairaut theory, noted with a (*) in 
Table [3l which corresponds to a value for H of: H* = 
3.26715240 x 10"^. 



3.3. Error estimation 

We can estimate the error that the use of the Clairaut the- 
ory introduces into the AH results. Indeed, if we consider 
the MHB value as the reahstic H value (see Table [3]), the 
relative error made is: 



H 



o-ff = 



MHB 



H* 



Hmhb 



2 X 10" 



(19) 



We estimate that the error is about 0.2 %. So, comput- 
ing the variable part of H with the C20 data results in a 
maximum error of about: 



lAiJrcai - AH*\ ~ (2 X 10-'"^) X (6 X 10-9) 



~ 1.2 X 10 



(20) 



assuming that the maximum value for AH is of the order 
of 6 X 10"^. Then, regarding the values of the AH data and 
of their precision, we can consider this error as negligible. 

4. Time series of AC20 used in this study 

The geodetic data used are the time series (variable part) 
of the spherical harmonic coefficient C20 of the geopo- 
tential, obtained by the GRGS (Groupe de Recherche en 
Geodesie Spatiale, Toulouse) from the precise orbit deter- 
mination of several satellites (like LAGEOS, Starlette or 
CHAMP) from 1985 to 2002 (Biancale et al., 2002). The 
combination of these satellites allows the separation of the 
different zonal geopotential coefficients, more particularly 
of J2 and J4. This series includes (i) a model part for 
the atmospheric mass redistributions (Chao & Au, 1991; 
Gegout & Cazenave, 1993) and for the oceanic and solid 
Earth tides (McCarthy, 1996), and (ii) a residual part (see 
Fig. [2]) obtained as difference of the space measurements 
with respect to a model. These various changes in the 
Earth system are modelled as variations in the standard 
geopotential coefficient C20 and we note the different con- 
tributions AC2o„t„, ^C'20oci ^C'2o,„,t,d ^^^d AC2o,e. : re- 
spectively. 

4.1. AC20 residuals and its secular trend: observed 
part 

Earlier studies already took into account the effect of the 
secular variation of C20 on the precession of the equator. 




Fig. 2. Normalized AC20 residuals (top: raw residuals, 
bottom: filtered residuals, where the high frequency sig- 
nals have been removed) : non-modelled part of the AC20 
harmonic coefficient of the Earth gravity field. 



Such a secular variation is attributed to the post-glacial 
rebound of the Earth (Yoder et al., 1983), which reduces 
its flattening. Wilhams (1994) and Capitaine et al. (2003) 
considered a J2 rate value of —3 x IQ-^/cy. Using the nu- 
merical value of Table [21 for the first order contribution 
(/oi coseo) to the precession rate rg, which is directly pro- 
portional to J2, the contribution J2/J2 x /oi coseo of the 
J2 rate to the acceleration of precession (PipA/dt^ is about 
— 0.014"/cy^, giving rise to a — 0.007"/cy^ contribution to 
the t^ term in the expression of tpA- 

Since 1998, a change in the secular trend of the J2 
data has been reported (Cox & Chao, 2000). This change 
can be seen in the series of AC'20 residuals (see Fig. [2]) . An 
attempt to model this effect, with oceanic data, water cov- 
erage data and geophysical models, has been investigated 
by Dickey et al. (2002). Using the residuals AC'20 of the 
GRGS, we can estimate a secular trend for J2 = — v^ 6*20 
from 1985 to 1998 (see Fig. E]). We find a J2 rate of the 
order of: —2.5 (±0.2) x lO^^/cy, which gives a change of 
about — 0.006" /cy^ in the t^ term of the polynomial de- 
velopment of the precession angle V'A ■ 

As this secular trend is not the same in the total data 
span, we will also model the long term variations in the 
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Table 3. Comparison between different values of the coefficient C j [MRe^) and of the constant part for H: (1) lAG 
values (Groten, 1999) - (2) MHB values (Mathews et al., 2002) - (3) Constant part H** obtained from Eq. JH]) using 
the M, Rp and C lAG values - (4) Method of "Clairaut" (Sect. 3.2), assuming hydrostatic equilibrium. The third and 
fourth methods use a constant part for C20 of -4.841695 x 10"'' in Eq. © (i.e. J2 = 1.0826358 x 10"^)^ ^he sense of 
the computation is indicated by the arrows. 





(1) 
lAG (1999) 


(2) 
MHB 2000 




(3) 

Separate values for 

M, 7?e and C 




(4) 
Clairaut 
Theory 


C/ (M Re^) 

H 


0.330701 
±2 X 10^'^ 

it 

3.273763 x 10"^ 
±2 X 10~* 


0.330698 

it 

3.27379492 x lO-^* 


H 


0.330722** 

* = 3.27355562 x lO-^* 


H* 


0.331370* 
= 3.26715240 x 10"^ 



C20 residual series with a periodic signal. Such a long- 
period term in the J2 residual series may come from mis- 
modelled effects, particularly from the 18.6-yr solid Earth 
tides. We will make such an assumption and adjust for the 
period 1985-2002, a secular trend and a long-period term 
in the AC20 residual series (see Sect. 4.3). 

However, it should be noted that a secular trend for J2, 
of the order of —3 x 10~^/cy, is more consistent with long 
term studies of the Earth rotation variations by Morrison 
& Stephenson (1997), based upon eclipse data over two 
millennia (they found J2 = (—3.4 ±0.6) x 10^^/cy). 

4.2. AC20 geophysical data used: modelled part 

The geophysical models that have been previously sub- 
tracted from the C20 data (i.e. atmospheric, oceanic and 
solid Earth tides effects) must be added back to these data 
in exactly the same way they had been subtracted to re- 
construct the relevant geophysical contributions. 

For each contribution we give the associated potential 
U at the point (r, (p, A, t) (limited to the degree 2 and 
order 0) that we identify with the Earth gravitational 
potential. Hence, we obtain the AC'20 coefficient contri- 
bution of each geophysical source. 

• The atmospheric contribution is due to pressure 
changes in time, measured and given by the European 
Centre for Medium-range Weather Forecasts (ECMWF) 
(see Fig. SJ. The simple-layer atmospheric potential, lim- 
ited to degree 2 and order 0, can be expressed as: 



£/«,,. = 47rGi?, 



1 + ^2 

5g 



AG 



20 



ECMWF 



where G = 6.672 x 10" 



kg- 



(i)P2o(sin0)(21) 



is the gravita- 



tional constant, ^2 == —0.314166 is a Love number (Farrell, 
1972), g ~ 9.81 m s^^, and P2o{sm(j)) is the Legendre 
function of degree 2 and order 0. The G, 



20 



ECMWF 



(t) atmo- 
spheric coefficient, expressed in Pascals, comes from the 
spherical harmonic decomposition of the ECMWF atmo- 




J2 filtered residuals 
Linear regression 




Fig. 3. J2 GRGS residuals (top: raw residuals, bottom: 
filtered residuals, where the high frequency signals have 
been removed): estimation of the linear trend, from 1985 
to 1998. 

spheric pressure grids, every 6 hours, over continents (see 
Gegout & Cazenave (1993) or Chao & Au (1991)): 



AC- 



20 



ECMWF ^ 



(i)- / Ap(0,A,i) 



dS{22) 



where S" is a surface grid pressure around the Earth and 
Ap is the difference of pressure with a constant part pre- 
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fixed, at the point {(p, A). Hence, identifying Eq. (PT|) with 
the Earth gravitational potential, we obtain the atmo- 
spheric pressure contribution to the AC2a harmonic co- 
efficient: 



AC,. 



,W 



47r Re^ (1 + fca) 
5Mg 



^^20ecmwf (*) 



(23) 



• The contribution of the oceanic tides (see Fig. [51) is 
modelled in the lERS Conventions 1996. The Earth re- 
sponds to the dynamical effects of ocean tides, and the 
associated potential, limited to the degree 2 and order 0, 
is: 

C/oc - 47r G i?e Pro ^-^ ( —) Ao(sin0) a{t) (24) 



where we note a, depending on time, as: 

a = 5m ^",2,0 cos(0„(f)+x„)+5±2,o sin(0„(t)+x„)(25) 

n + 

The sum over n corresponds to the Doodson develop- 
ment whose associated arguments are 0n and Xn- The 
parameter p^, (= 1025 kg m~^) is the mean density of 
sea water. Furthermore, C„ 2 o ~ ^n 2 sin(e„ 2 0) ^^'^ 

-5^,2,0 = ^n,2,o cos(e^,2,o)' '^hcrc (7^2,0 and e^2,o are the 
normalized amplitude and phase of the harmonic model 
of the oceanic tides limited to degree 2 and order 0. 
Identifying Eq. ([M)) with the Earth gravitational potential 
gives the oceanic tide contribution to the AC'20 harmonic 
coefficient: 



AC20oc(t) 



4^ j?e' (1 + fc^) PW 

5 M 



a{t) 



(26) 





1990 1995 

Years 



Fig. 4. Normalized atmospheric AC20 (top: raw data, bot- 
tom: filtered data, where the high frequency signals have 
been removed): atmospheric modelled part of the AC20 
harmonic coefficient of the Earth gravity field, obtained 
with ECMWF pressure data. 



• The solid Earth tide contribution (see Fig. [5]) is due 
to the gravitational effect of the Moon and the Sun on the 
Earth (lERS Conventions 1996). This force derives from a 
potential, developed in spherical harmonics, which limited 
to degree 2 and order is: 



B - 

UsoWid = G M — ^ P2o(sin0) C20 



Moon-I-Sun 



(t) (27) 



where 



C: 



20Moon+Sun 



{t) = 



, 3 sun 

fc20 ii-e 



5 M 



^ (^F2o(sin0,)j(28) 

'P—moon ^ P ^ 



where fc2o = 0.3019 is the nominal degree Love number 
for degree 2 and order 0, nip the mass of the body p, 
and Tp the geocentric distance and (j)p the geocentric lati- 
tude at each moment of the body p. The Love number de- 
pends on the tidal frequencies acting on the Earth. Hence, 
the contribution to AC20 from the long period tidal con- 
stituents of various frequencies v must be corrected (see 
lERS Conventions 1996). Eq. ^ corrected for the fre- 
quency dependence of the Love number, can be identi- 
fied with the Earth gravitational potential. We obtain the 



Earth solid tide contribution to the AC20 harmonic coef- 
ficient: 

^^203,iti<i(*) = ^20Moon+Su„ + "/^e9Mency correction" (29) 

This contribution comprises a constant part in the AC'20 
solid Earth tide, which is called "permanent tide". We 
have estimated it and obtained: —4.215114 x 10~^ (the 
lERS Conventions value is —4.201 x 10"^). We must re- 
move it from our AC20 data coming from solid Earth tides. 

• Finally, we must consider a series including all the 
effects described before. Hence, we add them back to the 
residuals (Fig. [2]), interpolating and filtering the data. 
Then we obtain the total series (Fig. [7]). 



4.3. Adjustments in AH data 

Eq. (jlip allows us to transform the geodetic AC20 tempo- 
ral variations into the dynamical flattening variations AH. 
They can then be introduced into the precession equations 
dl]), replacing H with {H + AH) (Eq. Q, Eq. Q, Eq. Q 
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Filtered and interpolated oeean tides C20 




Fig. 5. Normalized oceanic AC20 (top: raw data, bot- 
tom: filtered data, where the high frequency signals have 
been removed): oceanic-tide-modelled part of the AC20 
harmonic coefficient of the Earth gravity field; lERS 
Conventions 1996. 



Fig. 6. Normalized solid tides AC20 (top: raw data, 
bottom: filtered data, where the high frequency signals 
have been removed) : solid-Earth-tide- modelled part of the 
AC20 harmonic coefficient of the Earth gravity field; lERS 
Conventions 1996. 



and Eq. ([7])) and using the process already described in 
Sect. 2.3. 

It is generally considered that VLBI observations of 
the Earth's orientation in space are not sensitive to the 
atmospheric and oceanic contributions to the variations 
in C20 (de Viron, 2004). However the amplitudes of these 
effects have been evaluated in Table [TT] for further dis- 
cussion and in any case we can notice that they have a 
negligible effect on precession. 

The analytical and semi-analytical approach to solv- 
ing the precession-nutation equations provides polynomial 
developments of the tpA and to a quantities. The AH data 
are then considered as a linear expression plus Fourier 
terms with periods derived from a spectral analysis (18.6- 
yr, 9.3-yr, annual and semi-annual terms) (see Tables|4l[5l 
[6]and[7]). Note that the phase angles used for adjusting the 
AH periodic terms are those of the corresponding nuta- 
tion terms. This implies changes in the development of the 
equatorial precession angles {^a, ^a)i which we describe 
in the next section. 

For the residual contribution of AH, we will consider 
(i) an adjustment of a secular trend over the interval from 



1985 to 1998 (see Tableland Eq. dSO])), and (ii) an ad- 
justment of a secular trend plus a 18.6-yr periodic term 
(see Tableland Eq. pTjl ). both added to the seasonal 
terms. The fit (i) of the secular trend gives: 

H ~ -7.4 X lO^Vcy ^ J2 ~ -2.5 x lO'Vcy (30) 

and the fit of model (ii) gives: 

AH = (74 X 10~") X t + (20.9 X 10"") X sin(wt) 

+ (32.5 X 10"") X cos(wi) (31) 

with Lo = 27r/0.186. 

We must recall that these adjustments have been made 
together with the fit of annual and semi-annual terms. In 
contrast, the higher frequency terms appearing in the AH 
data have been filtered and we therefore did not take into 
account other contributions, as for example the diurnal 
effects of the geophysical contributions in AH. 

5. Effects of the AH contributions on the 
precession angles 

On the basis of the models fitted to the time series of AH 
in the previous section, obtained with geodetic AC20 se- 
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Table 5. Adjustment of periodic terms in the AiJ contributions, for the data span 1985-2002 for various AH geo- 
physical sources (atmospheric A7f atm , oceanic tides A7f oc and soHd earth tides /S.H soitid , as well as the residuals 

Ai?,es.)) 



Units are in 10 ^ 


Orad. 














period 


AH 


res. 


A-ff atm. 


AH 


oc. 


AH 


soitid. 


(in years) 


sin 


COS 


Sin COS 


sm 


COS 


sm 


COS 


1 


2.17 


-4.02 


0.96 -1.66 


-3.92 


1.41 


-4.64 


0.89 


0.5 


-0.43 


3.71 




0.76 


1.56 


-0.34 


28.04 


18.6 


2.09 


3.25 








-0.46 


-26.29 


9.3 


-0.17 


-1.65 








-0.01 


0.29 





Fig. 7. Normalized total AC20: top is the total series in- 
cluding atmospheric, oceanic tides and solid earth tides 
effects and the residuals; bottom is the total series with- 
out the sohd earth tides effect. 



ries, we investigate the influence of these geodetic data on 
the precession angle developments. First, we evaluate the 
effect of the secular trend considered in the AC20 residual 
series. Second, we report on the influence of each geophys- 
ical contribution, on the influence of the residuals and on 
that of the total contribution. Finally wc focus on the pe- 
riodic effects resulting from the various AH contributions. 



Table 4. Summary of the constant parts for H and C20 
(Constant part + Permanent tide) used in this study. 



Hmkb 




3.2737942 x 10"^ 


C20 




-4.841695 X 10"* 


J2 




1.0826358 X 10~^ 


H* 




3.2671524 x 10"^ 


H* with geophysical 


3.2671521 X 10"^ 


constant 


parts 




H** 




3.2735556 x lO^^* 



Table 6. Speciflc adjustment of the AH residual series 
(AiJrcs.), from 1985 to 1998. The secular trend is consid- 
ered as in Eq. ([30)l - Units are in 10^^" rad. 



period 
(in years) 



AH res. 

sin cos 



1 
0.5 



2.57 
-0.50 



Table 7. Adjustment of the total series of AH {AH tot.), 



from 1985 to 2002 - Units are in 10' 



rad. 



period 
(in years) 



AH tot. 
sin cos 



1 


-5.39 


-3.39 


0.5 


0.07 


33.67 


18.6 


0.92 


-23.11 


9.3 


-0.08 


-0.50 


influence 







5.1. Jo 



We have already mentioned that the J2 influence was 
taken into account in previous precession solutions 
(Wilhams 1994, Capitaine et al. 2003) (see Sect. 4.1). But 
depending on the value adopted, the polynomial devel- 
opment of the ipA precession angle is different. Indeed, 
if we take J2 = —2.5 x 10^^^/cy like in our study, or 
ij = -3 x 10""/cy like in Capitaine et al. (2003), the 
contribution in tJja varies by about 1.5 mas/cy^ (see Table 
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[5]). So we must carefully take into account this J2 rate. 
Furthermore, (i) we already noticed that such a secular 
trend has been recently discussed because of the change 
in this trend in 1998 (see Fig. [2) and (ii) the uncertainty in 
this secular trend, derived from space measurements of J2 , 
is significant. Therefore we can conclude that until there 
is a better determination of the J2 rate, the accuracy of 
the precession expression is limited to about 1.5 mas/cy^. 

5.2. Precession 

First, we can compare the polynomial part of our solu- 
tion Geod04 for the precession angles, based on the con- 
stant part Hmhb of H and on its variable part provided 
by expression ([?T|) . with previous precession expressions 
(IAU2000 and P03) (see Tabled. The differences larger 
than one /las concern the ipA precession angle and more 
particularly its t^ and t^ terms. The differences (of 7 mas 
and 2 /las, respectively) with respect to P03 are due to 
considering or not considering the J2 effect. Actually, P03 
includes a J2 secular trend, whereas Geod04 includes in- 
stead a 18.6-yr periodic term (see (3) in Table [9] or (2) in 
Table [TU]). Comparing Geod04 with the lAU 2000 preces- 
sion (which does not consider the J2 rate) shows differ- 
ences of 0.6 mas and 5 /ias in the t^ and t^ terms, respec- 
tively. This results from the improved dynamical consis- 
tency of the Geod04 solution (based on the P03 precession 
equations) with respect to lAU 2000. Note that such re- 
sults regarding the t^ and t^ terms will not be affected if 
changes of the order of 1 mas/cy in the precession rate 
would occur in an updated P03 solution. 

Second, we can evaluate the differences introduced in 
the ipA (and loa) polynomial development by the use of a 
constant part for H determined with the geodetic J2 (as 
used in GeodO^-H* and GeodO^-H**) instead of the Hastro 
determined by VLBI and used in Geod04. Table [5] shows 
that the differences are very large, but it should be noted 
that using J2 for deriving H suffers from the too large 
errors introduced by the mismodelled C/MRe^. 

5.3. Periodic contributions 

On the basis of the adjustments made in Sect. 4.3 for the 
different AH contributions, we estimate here the periodic 
effects appearing in the expressions of the precession an- 
gles. We can focus on the Fourier terms in the ijja preces- 
sion angle, which are the most sensitive to the AH effects. 
The corresponding results are presented in Table [TTl 

• First, we note that the major effect is due to the 
18.6-yr periodic term in the solid Earth tides (contribution 
number (3) of Table [TT|) : about —2 /j,as and 120 /ias in 
cosine and sine, respectively. The tidal annual and semi- 
annual effects are negligible as well as the atmospheric and 
oceanic effects (contributions number (4) and (5) of Table 

EI). 

• Second, note that the AH variation strictly limited 
to its residual part, not modelled into the geodetic orbit 



restitution, introduces negligible Fourier terms into the ipA 
development. But we can note that the way the long term 
effect is considered in such data (i.e. either with a secular 
trend term (contribution number (1) of Table [TT|) or a 
18.6-yr periodic term (contribution number (2) of Table 
llip ) is important. Modelling the long term variation in 
the geodetic residuals over the total data span as a 18.6- 
yr variation induces a term with an amplitude of 15 /zas in 
the ipA development. But at present the AC20 data span 
is not long enough to allow us to discriminate between the 
two models. 

Finally, we can conclude that the geodetic determina- 
tion of the total variable C20 (contribution number (6) of 
Table fTTj) introduces Fourier terms into the ipA precession 
angle development, mainly a 18.6-yr periodic one, of the 
order of 4 /ias and 105 /^as in cosine and sine, respectively. 



6. Discussion 

This study was based on new considerations: the use of 
a geodetic determination of the variable geopotential to 
investigate its influence on the developments of the pre- 
cession angles. The major effect on the precession is due 
to the J2 secular trend which implies an acceleration of 
the ^A precession angle. But for the moment, the avail- 
able time span for J2 satellite series is not as long as we 
need to determine a reliable J2 value. The J2 secular trend 
estimation based on our C20 residuals series from 1985 to 
1998 is: J2 = —2.5 x 10~^/cy. The accuracy of the preces- 
sion expression is limited to about 1.5 mas/cy^ due to the 
uncertainty in this J2 rate value. 

Then, we can notice that the main periodic effect is due 
to the 18.6-yr periodic term in AC2O due to solid Earth 
tides. But we must say that computing the AC20 with 
satellite positioning observations requires making some as- 
sumptions on the geophysical contributions to AC20, for 
instance from atmospheric pressure, and oceanic or solid 
Earth tides. Actually, models are used, but they are not 
perfect and we may have some errors. So the AC20 resid- 
uals obtained may be affected by these errors, which is 
why the total AC20 contributions (residuals observed -I- 
models assumed) constitute a better series to evaluate the 
effects on the precession angles. This introduces Fourier 
terms into the tpA development (4 /ias and 105 /ias in co- 
sine and sine respectively; see Table [TT]) that we should 
compare to the MHB2000-nutations. Indeed, the differ- 
ent terms of the total AH (or AC20) contributions have 
same periods as the (A";/;, Ae) nutations. This implies that 
there is some coupling between the observed AH effects 
and the nutations, which may not have been included in 
the MHB2000-nutations. 

In the future, we will be able to compare the J2 data 
with geophysical models and data, in order to have better 
ideas on the different contributions and on the secular 
trend. We will also be able to proceed to numerical study 
of this problem, and to implement a refined and more 
realistic Earth model. 
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Table 8. Influence of J2 on the polynomial development oiipA (more particularly on the t^ and t^ terms): (1) IAU2000 
(Mathews et al., 2002), (2) P03 (Capitaine et al., 2003) and (3) Same computation as in P03 but with other J2 values. 
The J2 secular trend estimation based on our C20 residuals series is: J2 = —2.5 x 10~^/cy. 





J2 


i^ 


t' 


(1) IAU2000 


None 


-1".07259 


-0". 001147 


(2) P03 


-3 X lO^Vcy 


-1". 079007 


-0". 001140 



Difi'erences wrt P03 



(3) 



0/cy 
-2 X lO^Vcy 
-2.3 X lO^Vcy 
-2.5 X lO^Vcy 



-7.000 mas 


2 /ias 


-2.871 mas 


1 /ias 


-1.954 mas 


1 pas 


-1.495 mas 


1 /ias 



Table 9. Polynomial part of the iJ^a and lua developments (units in arcseconds): comparison of (1) IAU2000 (Mathews 
et al. 2002) - (2) P03 (Capitaine et al. 2003) - (3) Differences of Geod04 (this study) with respect to P03, considering 
all the contributions for AH (Table El AH tot) - (4) Differences of Geod04 with respect to P03, obtained with a H 
constant part different from Hmhb, but not used in the following (see Table [3] for the H* and H** constant values). 



Angle 


Source 


i« 




t 


t' 


r' 




(1) 
(2) 


IAU2000 
P03 






5038". 478750 
5038". 481507 


-1".07259 
-1".079007 


-0".001147 
-0".001140 


Va 


(3) 
(4) 


Geod04 Hmhb 
Geod04 1 ^,. 








Differences wrt P03 




0" 

~ 10".23 
~ 0".37 


-7 mas 

-9.177 mas 
-7.079 mas 


2 /ias 

-3 /ias 
2 /ias 






(1) 
(2) 


IAU2000 
P03 


84381" 
84381" 


.448 
.406 


-0". 025240 
-0". 025754 


0".05127 
0".051262 


-0". 00772 
-0". 007725 




(3) 

(4) 


Geod04 Hmhb 
Geod04 1 ^1 






Differences wrt P03 in /las 




UJA 


















104 
3 




-35 
-1 
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Angle 
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t 




t' 






t' 
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-0" 


.001140 


V'A 


(2) Geod04 total contributions 

/o^ r. jn. -J 1 f 1985-1998 

(3) Geod04 residuals <^ 1935,2002 






Differences wrt P03 


in 


/las 















-7000 

-1495 
-7000 






2 

1 

2 
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(units in /xas). 







Periodic contribution for the t" 


term 


oi IpA 










fj,as 
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sin 




(1) 


Residuals 


Annual 


-1 


-1 






(1985-1998) 


Semi-annual 


- 


1 




(2) 


Residuals 


18.6-yr 


10 


-15 






(1985-2002) 


9.3-yr 


- 


4 








Annual 


-1 


-1 








Semi-annual 


- 


1 


AH 


(3) 


Solid tides 


18.6-yr 


-2 


120 


periodic 






9.3-yr 


- 


-1 


contributions 






Annual 


1 


- 








Semi-annual 


- 


3 




(4) 


Ocean tides 


Annual 


1 


- 








Semi-annual 


- 


- 




(5) 


Atmosphere 


Annual 


- 


- 




(6) 


Geod04 


18.6-yr 


4 


105 






total 


9.3-yr 


- 


1 






contributions 


Annual 


1 


-1 








Semi-annual 


- 
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